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Abundance stratification in Type la supernovae - III. 
The normal SN 2003du 
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ABSTRACT 

The clement abundance distributions in the ejecta of Type la supernova (SN) is stud- 
ied by modelling a time series of optical spectra of SN 2003du until ~ 1 year after the 
explosion. Since SN 2003du is a very normal Type la SN both photometrically and 
spectroscopically, the abundance distribution derived for it can be considered as rep- 
resentative of normal Type la SNe. We find that the innermost layers are dominated 
by stable Fe-group elements, with a total mass of ~ 0.2M Q , which are synthesized 
through electron capture. Above the core of stable elements there are thick 56 Ni-rich 
layers. The total mass of 56 Ni is 0.65M Q . The Si- and S-rich layers are located above 
the 56 Ni-rich layers. The dominant element in the outermost layers (M r > 1.1M , 
v > 13000 kms" 1 ) is O, with a small amount of Si. Little unburned C remains, with 
an upper limit of 0.016 M@. The element distributions in the ejecta are moderately 
mixed, but not fully mixed as seen in three-dimensional deflagration models. 

Key words: supernovae: general - supernovae: individual: SN 2003du - nucleosyn- 
thesis, abundances 



1 INTRODUCTION 

Type la supernovae (SNe la) are thought to be a thermonu- 
clear explosio n of a C+O white dw a rf (WD) in a close binary 
syste m (see iNomoto et al. | 1 1994 iHillebrandt fc Niemeveri 
2000). Since the ignition is triggered when the WD reaches 
a mass close to the Chandrasekhar limit, the properties of 
SNe la are expected to be homogeneous. SNe la are believe 
to synthesize a large amount of 56 Ni and other Fe-group el- 
ements during the explosion, and thus they are one of the 
primary supply sources of Fe-group elements in the Universe. 
In addition, the radioactive decay of 56 Ni results in an ex- 
treme luminosity of SNe la. This large luminosity, and their 
relative homogeneity makes it possible to use them as a stan- 
dardisable candles with whi ch distances can be measured in 
the distant Universe (e.g.. IPhillipsI Il993l : iRiess et alj 1 19961 
ll998l : |Perlmutter et alj|l999h . 

However, the explosion mechanism of SNe la is still 
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unclear. Although it is widely accepted that the explo- 
sion s tarts as a deflagrat ion, i.e., a subsonic burning 
flame (|Nomoto et al.l Il976l h the subsequent evolution of 
the flame is under debate: th e flame might rema in sub- 
sonic (deflagration model, e.g.. INomoto et al]|l984l h or ex- 
periences a transition to the supe r sonic regime (delayed 
detonation model, e.g., iKhokhlovl 1 199lh . These scenar- 
ios have been studied in detail with multi-dimensional 



simulations (e.g.. |Reinecke et al. 


12002): iPlewa et all 


2004 Ropke & Hillebrandtl 2005: 


Gamezo et al. 2005; 


Ropke & Niemevei 2007: Jordan et all 20081: Kasen et all 


2009|). 



Different explosion models predict different element 
abundance distributions and nucleosynthesis yields. For ex- 
ample, deflagration models tend to leave more unburned 
elements than delayed detonation models, where the WD 
is burned almost entirely. Three-dimensional deflagration 
models predict a mixed abundance distributio n, which is not 
seen in one-dimensional deflagration models (|Nomoto et al.l 
1 19841 ) or in one- and multi-dimensional delay ed detonation 
models (|Khokhlovllr99ll : iGamezo et al J 120051 ) . Thus, in or- 
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der to place constrains on the explosion mechanism, it is 
important that the abundance distribution and the nucle- 
osynthesis yields are derived from observed SNe. 

Modelling a time series of optical spectra of extragalac- 
tic SNe is a powerful method to extract the abundance dis- 
tribution in SNe la. At early epochs (photospheric phase, 
iexp < 40 days, hereafter £ oxp denotes the time since the 
explosion), the optical spectra of SNe la exhibit a pseudo- 
continuum emission and P Cygni profiles superposed onto it. 
Since the absorption features are formed by line scattering 
in the outer, optically thin layers, they carry information on 
the abundances in the outer layers. In addition, since the 
SN ejecta are expanding, the line-formation region moves 
inwards with time in the mass and velocity coordinate. Fol- 
lowing the time evolution of the spectra, the abundances in 
different layers can thus be derived. 

At late enough epochs (nebular phase, t cxp > 150 — 200 
days), a SN shows an emission- line spectrum. Most lines 
are forbidden lines of heavy elements. Since these lines are 
formed in the innermost layers of the ejecta, we can derive 
the abundances in the inner layers, which cannot be studied 
with the photospheric-phase spe ctra. 

Using the optical spectra, IStehle et al.l (|2005l . Paper 

I) derived the abundance stratification in SN 2002bo. SN 
2002bo is a normally luminous Type la SN (with a B-band 
light curve (LC) decline rate Amis = 1.13 mag), but has 
higher line velocities, i.e., larger Dopple r shifts at the absorp - 
tion minima, than canonical SNe la |Benetti et al.|[2004h . 
Paper I found that a large abundance of 56 Ni and interme- 
diate mass elements (IMEs) such as Si and S is required in 
the outer layers, and pointed o ut that mixed - out 56 Ni can 
explain the fast rise of the LC. IMazzali et al.l (|2008l . Paper 

II) studied SN 2004eo. SN 2004eo has a transitional lumi- 
nosity between normal and subluminous SNe la and has a 
relati vely rapid decline rate of the light curve (Amis = 1-46 
mag. IPastorello et al.l 120071 ). A relatively small mass of Fe- 
group elements was obtained for SN 2004eo in Paper II, sug- 
gesting that this is responsible for the rapid decline of the 
LC. 

In this pa per, we study the abundance distribution 
in SN 2003du (ICerardv et all I200I lAnupama et al.l 120051 ; 



IStanishev et al.ll2007l ; iMarion et alj|2009h . SN 2003du is a 
very normal Type la SN both photometrically and spec- 
troscopically, and it can therefore be considered as a rep- 
resentative of normal SNe la. SN 2003du is classified as 
low velocity g radien t (LVG) according to the scheme of 
iBenetti et al] (I2005TI, or as co r e-norm al according to the 
scheme of branch et all (|2006l . 120071 '). This SN is an ex- 
tremely well-studied object at optical and near infrared 
(NIR) wavelengths, with an extremely good time sampling, 
and one of the rare cas es in which NIR spec t ra at late epochs 
have been obtained (|Hoflich et al.l |2004| ; iMotohara et al.l 
12006ft . 

This paper is organized as follows. In Section [2] the 
methods of analysis are described. The spectra at the pho- 
tospheric and nebular phases are modelled in Sections [3] and 
[4] respectively. In Section [5] the derived abundance distribu- 
tion is tested against the bolometric LC. The derived abun- 
dance distribution and the integrated yields are shown in 
Section [6] In Section [7] we give conclusions. 

Throughout the paper, th e distance to SN 200 3du is 
assumed to be ^ = 32.79 mag (jStanishev et all 120071 ). The 



reddening in our Galaxy and in the host g alaxy are E(B — 
V) = 0.01 mag and mag, respectively (|Stanishev et al.l 
l2007h . We adopt a bolometric rise t ime of 19 days, which 
was derived from the LC synthesis (|Stanishev et al.ll2007l ) 
and is consistent with the B-band rise time of 20.9 days. 

The mass and mass fraction of radioactive 56 Ni 
[M( 56 Ni) and A"( 56 Ni), respectively] are expressed by the 
value at the explosion (t cxp = 0). M(Fe) [A(Fe)] and M(Ni) 
[X(Ni)] only represent the mass [mass fraction] of stable Fe 
and Ni, respectively. 



2 METHODS OF ANALYSIS 
2.1 Photospheric Phase 

The photospheric-phase spectra are studied using the 
Monte Carlo s pectru m synthesis code developed by 
IMazzali fe Lucvl (|l993T ). Assumptions made in the code 
are described in Ap p endix |A"T1 F o r more detail s, see 
IMazzali fe Lucvl (|l993l ); iLucvl (|l999l ); iMazzalil (|2000l ). The 
code has been used for the modelling of various types of 
SNe (e.g., IMazzali et aT]|l992l . ll993l . l2000l ). 

The code calculates a synthetic spectrum for given in- 
puts: (1) a radial density profile, (2) a position of the photo- 
sphere in velocity (i.e., photospheric velocity w p h)El (3) an 
emergent luminosity L and (4) element abundances above 
the photosphere. F or the density struct ure, we use the defla- 
gration model W7 (|Nomoto et al.lll984h . The validity of this 
choice is discussed in Section [5] The other parameters are 
optimized so as to give the best fit of the observed spectra. 

The luminosity L is most strongly constrained by the 
observed absolute flux. The photospheric velocity v p h is esti- 
mated primarily by using the overall colour of the observed 
spectra via the relation L — 47™ ph fc Xp o"T ( ? ff , where T e g and a 
are the effective temperature and the Stefan-Boltzmann con- 
stant, respectively. The estimated velocity is further checked 
based on the reproduction of the observed line velocities (i.e., 
Doppler shift at the absorption minimum of P Cygni pro- 
files, see Fig. [TJ . 

With the estimated L and w p h, the abundance ratios 
above the photosphere are optimized to reproduce the ob- 
served spectra. We start modelling from earlier spectra for 
abundance stratification because the photosphere recedes 
with time in the mass and velocity coordinates. Suppose 
that the photosphere recedes from vi to V2 (vi > V2) in the 
time interval from t\ to t2 (ii < tz). First, the abundances 
at v > Vi are optimized by fitting the observed spectrum 
at f cxp = ti. At the later phase (t cxp = *2, v P h = V2), the 
inner part begins to make a contribution to the spectrum. 
However, even at t exp = i2, the outer layers at i) > i>i still 
contribute to the spectrum. Thus, we keep the abundance 
ratios at v > Vi and optimize the abundances at V2 < v < vi 
by fitting the observed spectrum at i cxp = £2- Fitting the 
time series of spectra, we get an abundance distribution as 
a function of velocity. 



Since the SN ejecta are in homologous expansion (v = r/i) 
at the time of observations, the velocity can be used as a radial 
coordinate. As the velocity is unchanged with time, it is more 
convenient than the radius itself. 
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Table 1. Parameters of the synthetic spectra at photospheric 
phases and the calculated effective temperature and blackbody 
temperature 
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2.2 Nebular Phase 



A nebular-phase spectrum is modelled under different as- 
sumptions from those made for the photospheric-phase spec- 
tra. The assumptions made in our anal ysis are describ ed in 
Appendix | A2 For mo re details, see e.g., lAxelrodl (|l980l ) and 
iMazzali et al. (|2001bh . Thi s code has been used for the mod- 
elling of supernovae (e.g., IMazzali et aT]|2001bt ). 

The inputs are only the radial density profile and the 
element abundance distribution. For the density structure, 
we use the W7 model as in the photospheric phase. For 
the abundance distribution, the results of the photospheric- 
phase modelling are used in the outermost layers (v > 15000 
kms" 1 ). At intermediate layers (8200 kms -1 < v < 15000 
kms -1 ), the abundance distribution is constrained both 
from photospheric-phase spectra and nebular-phase spectra. 
The abundances in the inner layers (v < 8200 kms -1 ) are 
constrained only from the nebular-phase spectra. The mod- 
elling is iterated until a consistent distribution is obtained. 

The total mass and distribution of 56 Ni are especially 
important for the modelling of nebular spectra because they 
determine the emergent luminosity and the energy deposi- 
tion in each layer. The total energy emitted by forbidden 
lines of Fe is sensitive to the mass of 56 Ni as well as to the 
total amount of Fe (i.e., the sum of stable Fe and the decay 
product of 56 Ni), which affects in particular the ionization of 
Fe. Thus, although 56 Ni is mostly Fe at the nebular epoch, 
stable Fe and 56 Ni products can be distinguished, albeit only 
indirectly. 



3 PHOTOSPHERIC-PHASE SPECTRA 

T welve photospher i c-phas e spectra of SN 2003du presented 
bv lStanishev et all (|2007| ') are modelled. The input param- 
eters (L and v p h) are summarized in Table [TJ The effec- 
tive temperature T c a is calculated by the equation L = 
47rUp h £exp0"r e 4 ff ■ The blackbody temperature Tb calculated 
with the code is also shown in Table [T] The blackbody tem- 
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perature is defined by Lb = *iTrv ph t cxp aT B , where the lu- 
minosity Lb in this equation is not the emergent luminos- 



Figure 1. Line velocities meas ured in the observed spectra (open 
symbols, IStanishev et al.ll200?f ) compared with the photospheric 
velocities derived from the spectral modelling (filled red circles). 
The photospheric velocity follows the Si 111 line velocity until 
~ 20 days after the explosion (around maximum). 



ity but the luminosity emitted from the photosphere, in- 
cluding photons scatt ered backward into the photosphere 
(IMazzali fe Lucvlll993h . Thus, Tb is always higher than the 
effective temperature (T e g). 

In Fig. [T] the position of the photosphere is shown 
by red circles. The observed line velocities of Si II A6355, 
Si ill A4560, and S u A5640 are also shown for comparison 
(|Stanishev et al.ll2007h . The time evolution of the estimated 
photospheric velocity is similar to that of the Si in line ve- 
locities, which is the slowest among the three lines. This 
indicates that the slowest line velocity is the best observa- 
tional indicator of the photospheric velocity. Note that the 
estimated photosphere after maximum (t cxp > 21 days) is 
not very certain because the assumptions of the code become 
less reliable at later epochs (see Appendix I A 10 , 

In the following sections, the results of modelling at pre- 
maximum epochs (from —12.8 to —5.8 days since B maxi- 
mum), maximum epochs (from —4.0 to +1.2 days) and post- 
maximum epochs (from +4.3 to +17.2 days) are presented. 



3.1 Pre-Maximum Epochs 

In Fig. [2] the observed spectra at epochs from —12.8 to —5.8 
days since B maximum (black) are compared to the syn- 
thetic spectra (red). The observed spectra at these epochs 
show a relatively well defined pseudo-continuum and P 
Cygni profiles of IMEs and Fe-group elements. At these 
epochs, the photospheric velocity decreases from 11500 to 
9400 kms" 1 . 

Fe-group elements: The presence of the Fe, Co and Ni 
lines clearly indicates that some amount of these Fe-group 
elements are located in the outer layers. Especially, since 
only a small amount of 56 Ni has decayed into Fe at such 
early epochs, the presence of Fe line s requires stable Fe in 
the outer layers jTanaka et al.|[2008h . The mass fraction of 
stable Fe is found to be X(Fe) ~ 0.003 - 0.005. 
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4000 6000 8000 10000 4000 6000 8000 10000 

Rest wavelength (A) Rest wavelength (A) 

Figure 2. The observed spectra of SN 2003du at the prc-maximum epochs (black) compared with the synthetic spectra (red). Panels 
(a) - (d) show the spectra at —12.8, —10.9, —7.8 and —5.8 days from B maximum, respectively. Line identifications shown in the figures 
are made only for the major contributions. The synthetic spectra are reddened with E(B — V) = 0.01 mag. 



Calcium: In the spectrum at —10.9 days from B max- 
imum, the Ca n lines near 4000 A (Ca II H&K) and 
8000 A (Ca II IR triplet) show an extremely high ve- 
locity (v ~ 25000 kms -1 ). The origin of this line has 
been discussed, including the interaction with the cir- 
cumstellar matter (CSM) or high velocity blobs (e.g. 



Hatano et al. 1999: Thomas et al. 2004; Gcrardv et al. 2004: 



Mazzali et all l2005bjlal: iQuimbv et all l200fj: ITanaka et all 



200d : lAItavilla et alj|2007i ; iGaravini et alj|2007t ). In this paT 
per, we do not consider the density enhancement with re- 
spect to the W7 model. Keeping the W7 density struc- 
ture, a mass fraction as high as Jf(Ca) > 0.3 is required 
at v > 22000 km s -1 . 

The density at the high velocity layers at v > 22000 
kms -1 in delayed detonation models is higher than in 
the W7 model by a factor of ~ 10 (s ee e.g., Fig. 11 of 
iBaron et af1l200uf ). ITanaka et all l|2008l ) studied the influ- 
ence of the density structure at the outermost layers. They 
found that the required mass fraction of Ca is decreased by a 



factor of ~ 10 in the delayed detonation models, but it is still 
higher than solar abundance (see their Table 3). The den- 
sity enhanceme nt may also possibly be caused by interacti on 
with the CSM (|Gerardv et alj|2004 iMazzali et al.ll2005bl ). 

Magnesium, Silicon, and Sulphur. Since the strongest 
Mg feature at 4200 A is contaminated by Fe lines, the mass 
fraction of Mg is not as well constrained as those of Si, S, and 
Ca. If we assume that the Fe abundance is constrained from 
the features around 4700 A, and we attribute any differences 
to Mg, the mass fraction of Mg is X(Mg) ~ 0.01 - 0.1 at 
v > 9400 kms" 1 . 

SN 2003du has relatively narrow Si II lines. To repro- 
duce the profile of the Si II A6355 line, the mass fraction of 
Si is required to be A" (Si) ~ 0.3. The dominant element in 
the layers at v > 13000 kms - is not Si, but more likely O. 
The distribution of Si is discussed in Appendix [B] 

The mass fraction of S is X(S) = 0.1 — 0.2 around v = 
10000 kms -1 . Although the synthetic spectra for the two 
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Figure 3. Same as Fig. [2] but at the maximum epochs. Panels (a) - (d) show the spectra at —4.0, —1.9, and +1.2 days from B 
maximum, respectively. 



earliest epochs show stronger lines than the observations, 
the later spectra are nicely reproduced. 

Carbon and Oxygen: There is no clear C line in the 
observed spectra except for the possible C II A6578 line 
at the emission peak of the strongest Si II A6355 line 
jStanishev et al.l 120071 ) . We derive an upper limit for the 
C mass as M(C) < 0.016 M© at v > 10500 kms" 1 (details 
are presented in Appendix [Cj. The small mass fraction and 
mass of C are consistent with the results of previous works 
(IMarion et al 1120061 ; iThomas et "ai]|2007l ; iTanaka et aLll200Sl ; 
iMarion et al"[|2009l ). 



The synthetic spectra show a clear O I A7774 line 
around 7500 A while the presence of such a strong O line is 
not that clear in the observed spectra. Although the abun- 
dance of O is not well constrained, considering the lower 
abundance of C, Mg, Si and S (see above), we conclude that 
the dominant element in the outer layers is O. The mass 
fraction of O at O 9400 kms" 1 is X(0) = 0.05 - 0.85, 
increasing towards the outer layers. 



3.2 Maximum Epochs 

In Fig. |31 the observed spectra at the epochs from —4 to 
+ 1.2 days since B maximum (black) are compared to the 
synthetic spectra (red). The element features present in the 
observed spectra are similar to those in the pre-maximum 
spectra. At these epochs, the photospheric velocity decreases 
from 9050 to 8200 kms^ 1 . The flux of the synthetic spectra 
at the line-free region at A > 6500 A tends to be overes- 
timated because of our assumption of blackbody emission 
(see Appendix (AT) . 

Fe-group elements: Around maximum, lines of Fe-group 
elements are present, as in the pre-maximum epoch. A larger 
amount of 56 Ni is required compared with the pre-maximum 
epochs to block the emission in the near-UV more effectively. 
The mass fraction of 56 Ni is X( 56 Ni) = 0.5 - 0.7 at v = 
8200 - 9050 kms -1 . 

At maximum epochs, the presence of Fe lines does not 
directly indicate the presence of stable Fe because of the 
high abundance of the decay product of 56 Ni. At B max- 
imum (t cxp = 20.9 days), about 10 % of 56 Ni has already 
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Figure 4. Same as Fig. [2] but at the post- maximum epochs. Panels (a) - (d) show the spectra at +4.3, +7.2, +10.0 and +17.2 days 
from B maximum, respectively. 



decayed into Fe, and thus an Fe abundance of ~ 0.05 — 0.07 
is the product of the decay of 56 Ni. This is much larger than 
the abundance of stable Fe required for the pre-maximum 
spectra, and thus we cannot distinguish the contribution of 
stable Fe. 

Silicon and Sulphur: The Si II A6355 line becomes nar- 
rower. The Si distribution at the outer layers derived from 
the line profile in the pre-maximum spectra (Appendix [B| 
nicely reproduce the maximum spectra. In the spectrum at 
maximum, the line ratio of two Si n lines (A5972 and 6355, 
ft(Si n) in iNugent et ail Il995h is also nicely reproduced, 
suggesting the estimates of the temperature and of th e ion- 
ization are reasonable (Table IT1 lHachinger et al1l200of ). 

The S II lines around 5500 A are also reproduced 
nicely. Although the photospheric velocity at this epoch is 
w p h= 8200 — 9050 kms" 1 , the S II lines, as well as the Si II 
lines, are mainly formed at v ~ 10000 kms -1 (see Fig. [T]). 
The "detachment" of the lines is caused by the ionization 
effect, e.g., S ill is dominant near the photosphere while the 



frac tion of S H increase s towards higher velocities (see Fig. 
9 of lTanaka et alj|2008f ). 



3.3 Post-Maximum Epochs 

In Fig. 31 the observed spectra at epochs from +4.3 to +17.2 
days since B maximum are compared to the synthetic spec- 
tra. At these epochs, the contribution from Fe-group ele- 
ments is larger than at earlier epochs. The assumption of 
blackbody emission becomes unreliable, and results in an 
increasing overestimation of the flux at A > 6500 A. Al- 
though the photospheric velocity decreases from 7700 to 
4000 kms -1 in our models (Table[T}, these values are there- 
fore not strongly constrained. 

Since the assumption that there is no energy deposition 
above the photosphere is also unlikely to hold at this epoch, 
we use the abundance distribution derived from the mod- 
elling of the nebular spectrum (Section [4}, rather than op- 
timize the abundances using the spectra at post-maximum 
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epochs. Even with this treatment, the synthetic spectra give 
reasonable agreement with the observed spectra. 

Fe- group elements: The 56 Ni distribution derived from 
the nebular spectrum reproduces the Ni, Co and Fe lines 
present in the observed spectra. The Ni II, Co n and Fe n 
lines are responsible for the strong blocking of the UV flux. 
The layers below 8200 kms -1 down to 4000 kms -1 (the 
photospheric velocity of the last spectrum, Fig[4ji) are rich 
in 56 Ni, which has a mass fraction of X( 56 Ni)=0.7 - 0.8. 
Note that the Fe lines at 4500 - 5000 A in the synthetic 
spectrum at t = +17.2 days are too strong. The discrepancy 
is also seen in the model for SN 2004eo (Paper I), and it 
may result from our assumption of no energy deposition in 
the atmosphere (see above). 

Sodium, Silicon, Sulphur, and Calcium: In the last spec- 
trum (Fig. 2Ji) , there is a strong absorption at 5600 A in the 
observed spectrum, which is expected to be the Na I D line. 
With our code, however, this line is not reproduced even if 
we set X(Na) =0.1. This is already too high compared with 
the prediction by the explosi on models: the yield of Na is 
only ~ 10~ 5 - 10" 4 M Q (e.g., Ilwamoto et~al]ll999l ). This is 
the result of the overestimate of Na ionization by our Monte 
Carlo code, which is a problem also for other NLTE codes 
l|Mazzali et al.|[l997l . and references therein), and it may be 
caused by inaccuracies in the atomic data. Thus, we do not 
discuss total mass of Na in the following section (Section[6]). 

The strong Si II and S n lines are reproduced well in the 
spectra at +4.3 and +7.2 days from B maximum (Fig. [3^ 
and[4p). These lines are formed in layers detached from the 
photosphere (see Fig. [TJ because of the high mass fractions 
of Si and S, and the large fraction of Si n and S n in layers 
above the photosphere. 

The Ca n H&K and IR triplet lines in the synthetic 
spectra show too strong high velocity component. They are 
caused by the high abundance at v > 15000 kms -1 derived 
from the earliest spectra. We believe that the constraints 
from the earlier spectra are more reliable because the valid- 
ity of the assumptions in the atmosphere is more certain at 
those epochs. 




4000 6000 8000 

Rest wavelength (A) 

Figure 5. The nebular spectrum of SN 2003d u (+376.9 days from 
B maximum, black line. lStanishev et all2007h compared with the 
synthetic spectra (color lines). The red solid line shows the best 
fit model while the blue dashed line shows the model without 
central region devoid of 56 Ni. 

If a central region devoid of 56 Ni is absent (blue dashed 
line), both Fe features look very sharp, unlike the observa- 
tions. In order to reduce the flux, the 56 Ni distribution must 
be limited to a smaller region, and the overall higher density 
of the inner region leads to a lower Fe ionization, so that the 
ratio of the two Fe lines in the blue is no longer correct: the 
Fe n-dominated line (near 5200 A) becomes too strong, as 
shown in the model in blue dashed line in Figure [5] On the 
other hand, the presence of stable Fe acts as a coolant, keep- 
ing the ionization degree to a level which yields reasonable 
line intensities (red line). 



4 NEBULAR-PHASE SPECTRA 

A spectrum at + 376.9 days since B maximum 
l|Stanishev et al.ll2007T ) is modelled to derive the abundances 
in the deepest layers. A 56 Ni distribution with a decreasing 
mass fraction toward higher velocities (v > 10000 kms -1 ) 
nicely fits the profile of the two strongest Fe lines (red line 
in Figure [5| • This is consistent with the constraints from 
the spectra at the pre-maximum and maximum epochs. 

The innermost layers are dominated by stable Fe. In 
particular, stable Fe and Ni dominate the abundance at ve- 
locities inside of about 3000 kms -1 , accounting for a total 
mass of ~ 0.271/0. Most of this is actually stable Fe. Lim- 
its to the stable Ni mass can be set by the weakness of the 
[Ni n] emission near 7400 A. As for stable Fe, this is required 
in the innermost regions in order to keep an ionization bal- 
ance between Fe II and Fe in which allows the shape of the 
two strong, Fe-dominated emissions at 4500 — 5000 A to be 
reproduced, as well as their ratio. The bluer emission, near 
4600 A, is in fact dominated by Fe in, while the redder one, 
near 5200 A, is dominated by Fe II. 



5 BOLOMETRIC LIGHT CURVE 

The derived ab undance distribution is tested against the 
bolometric LC (|Stanishev et alj|2007l ). SN 2003du is an ex- 
tremely well-observed SN, and the bolometric LC was con- 
structed wit h a small uncertainty, including the contribution 
of the NIR (jStanishev et al.ll2007r i. 

For the modellin g, we use a Monte Carlo LC code based 
on that developed in lCappellaro et al.l (|l997f h The code fol- 
lows the emission and propagation of 7-rays and positrons 
from 56 Ni and 56 Co decay. Both 7-rays and positrons are 
treated with an effective opacity (i.e., the positrons are 
not assumed to deposit in situ). The opacities adopted 
are 0.027 cm 2 gm -1 for the 7-rays and 7 cm 2 gm _1 for 
the positrons (|Axelrodl Il980l ) . For o ptical opacity , we use 
an a nalytic formula introduced by iMazzali et al.l (|2001al . 
l2007h . which includes the composition dependence (Fe- 
group elements and other el ements) and a time dependence 
iHoeflich fc Khokhlovlll996h . 

Fig. [S]shows the comparison between the observed bolo- 
metric LC (black squares) and the computed LC (blue line). 
The input luminosities for the photospheric-phase spectra 
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Days after the explosion 

Figu re 6. Bolomet ric LC of SN 2003du (black squares, 
IStanishev et al] 12007!) compared with the bolometric LC model 
(thin line) and the input bolometric luminosities in the spectral 
modelling (red circles). 

are also plotted (red circles). The model LC matches well 
with the observed LC from before maximum until > 1 year. 
This suggests that the derived abundance distribution, espe- 
cially that of 56 Ni, is reasonable. The bolometric LC derived 
from the synthetic spectra at epochs after maximum devi- 
ate progressively from the observational points, reflecting 
the increased production of spurious flux redwards of 6500 
A in the synthetic spectra (see Sect. 4). 



6 ABUNDANCE DISTRIBUTION AND 
INTEGRATED YIELDS 

6.1 Abundance Distribution 

The abundance distribution derived from the modelling of 
the photospheric- and nebular-phase spectra is shown in 
Figs. [7] and [8j These two figures show the mass fraction of 
each element as a function of mass and velocity, respectively. 
For compari son, the abundance d istribution of a deflagration 
model (W7. lNomoto et al.lll984l ) is also shown (Figs[7p and 
[Hp). Note the line with " 56 Ni" represents the mass fraction 
of radioactive 56 Ni at the time of the explosion (£ eX p = 0), 
and that "Fe" and "Ni" denotes only stable elements. 

At the innermost layers (M r < O.1M , v < 3000 
kms -1 ), stable Fe dominates (Section^. Th is is also sug- 
gested by the NIR spectru m at late phase |Hoflich et al.l 
12004 iMotohara et all 120061 ). The existence of stable ele- 
ments at the innermost layers is consistent with the ex- 
plosion model. These elements are synthesized via electron 
capture reactions at high densities, and thus are typically 
located in the innermost layers. 

Above the layer of stable elements, there is a 56 Ni - 
dominated layer (M r = 0.1 - O.9M , v = 3000 - 10000 
kms" 1 ). This is fairly consistent with the distribution of 
W7. However, in SN 2003du a moderate amount of 56 Ni is 
distributed also in the outer layers, which is important for 



the blocking of UV light in the early spectra. This was also 
true for SN2002bo (Paper I). At M r > O.9M (v > 10000 
kms -1 ), there is ~ 0.01M Q of 56 Ni. 

The layers dominated by Si and S are located above 
the 56 Ni -rich layers. Compared with W7, Si has a broader 
distribution in the mass/velocity coordinates. The moderate 
amount of Si in the 56 Ni -rich layers is required to explain 
the Si absorption line in the spectra after maximum. Even 
at the outer layers at v > 15000 kms -1 (M r > I.3M0), 
X(Si) = 0.01 - 0.1 is also required (Appendix[B|. This is in 
contrast to the W7 deflagration model, where C+O remain 
unburned at v > 15000 kms -1 . 

Above the Si and S-rich layers, the dominant element is 
O (M r > I.IMq, v > 13000 kms" 1 ). Although this is not 
directly constrained, a high abundance of Si would not be 
consistent with the observed line profile (Appendix |B]) and 
the mass fractions of C, Mg, S, and Ca are also not larger 
than O. 

The mass of unburned C is quite small, M(C) < 
O.O16A/ at v > 10500 km s -1 (Appendix [CJ. The small C 
mass indicates that the progenitor C+O WD was almost 
completely burned. In the outermost layers (M r > 1.3Mq, 
v > 15000 kms -1 ) only C-burning occurs, and elements 
heavier than O are not synthesized effectively. Note that 
there is no constraint on C above v = 18000 kms -1 , and 
a mass fraction as high as X (C) = 0.5 can be allowed, (see 
Appendix IC)) . 

It is also possible that the C abundance in the outer- 
most layers of the progenitor WD is already suppressed in 
the pre-SN stage. In the single-degenerate scenario, the WD 
grows by mass accretion from a companion star. In the case 
of weak He flashes, their nucleosynthesis products can ac- 
crete onto the WD. Since the products of the He flashes 
are no t always C and O but c ould also be heavier elements 
(e.g., lHashimoto et al1ll983l ). the outermost layers of the 
WD do not necessarily consist only of nearly equal mass 
fractions of C and O. The heavier elements synthesized in 
the pre-SN stage can also be responsible for the high velocity 
features of the Si an d Ca lines present in the earliest spectra 
(|Tanaka et alj|200gf ). 

The overall abundance distribution is not consistent 
with the almost fully mixed distri bution seen in three- 
dimensional deflagration models (|Travaglio et al.l |2004| ; 
iRopke fc Hillebrandtl |2005l ). But the boundaries between 
stable Fe/ 56 Ni, 56 Ni/IME, IME/O-rich layers are less clear 
than in the one-dimensional models. T his requires moderat e 
degree of mixing during the explosion (|Wooslev et al.ll2007l ) . 

The sensitivity of the results on the parameters has been 
studied in Paper II. In the photospheric phase, the synthetic 
spectra are sensitive to the velocity at the inner boundary 
and the element mass fractions. Variations in the mass frac- 
tions of up to 25 % are possible, while the synthetic spectra 
are more sensitive to variations of the velocity: a change of 
more than 15% can noticeably affect the fits (see Paper II). 
In the nebular phase, no given inner boundary is required 
and the constraints on the mass fractions are less uncertain. 

Except for the elements shown in Figures [7] and |5J 
there must be some amount of Ti and Cr in the ejecta. For 
photospheric-phase spectra, these elements are responsible 
for the features or the flux level around 3000 — 3500 A. The 
mass fraction of the sum of Ti and Cr is about 0.001 — 0.01 
(Papers I and II). However, since this wavelength range is 
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Figure 7. Abundance distribution in the mass coordinate derived from the spectral modelling (a) compared with that of W7 (b). " 56 Ni" 
represents the mass fraction of 56 Ni at t oxp = and "Fe" and "Ni" only represent mass fractions of stable Fe and stable Ni, respectively. 
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Figure 8. Same as Fig. [7] but in the velocity coordinate. 
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Atomic Number Z 

Figure 9. Abund ance ratios, i.e., [El ement /Fe] = log 1 n \M /M(Fe)} - log ln fX/X(Fe) |pi of SN 2003du (re d circles), W7 (black line, 
iNomoto et al.lll984t) , WDD1 (blue line. Ilwamoto et alj[l999h . and WDD2 (purple line. Ilwamoto et a.1.1 ^999^ ■ Note that, in this figure, 
the yields after the radioactive decay are shown, i.e., Fe includes both of stable Fe and the decay product of 56 Ni, and Ni includes only 
stable Ni. 



not covered in most spectra used in this paper, we refrain 
from discussing the distribution and total yields of these 
elements. 



6.2 Integrated Yields 

The ejected mass of each element in SN 2003du was cal- 
culated by integrating the abundance distribution assuming 
the W7 density structure. Results are summarized in Ta- 
ble [2] For comparison, the nucleosynthetic yields of the fol- 
lowing mode ls are also shown: on e-dimensional deflagration 
model (W7. INomoto et al.|[l983 ). one-dime nsional delayed 
deton ation models (WDD1 and WDD2 0, Ilwamoto et all 
19991). three-dimens ional deflagration model (b30_3d_768, 
Travaglio et alj|2004h . t wo-dimension al detonating failed de- 
flagration model (Y12, IPlewal [20071 ) . and tw o-dimensional 
off-cen ter delayed detonation model (O-DDT, iMaeda et al.l 
l2010bh . 

In Figure the derived yields are compared with those 
of the one-dimensional models in the abundance ratio to Fe, 
normalized by the solar abundances, i.e., [Element/Fe] = 
log 10 [M/M(Fe)] - log lo [X/X(Fe)] . Here, the yields after 
the radioactive decays are shown, i.e., Fe means the sum of 
stable Fe and the decay product of 56 Ni, and Ni includes 
only stable Ni. 

In SN 2003du, the ejected mass of neutron-rich elements 
is estimated to be ~ Q.2Mq. They are synthesized by elec- 
tron capture reactions and mostly distributed in the inner- 
most layers. The mass of 56 Ni is 0.65Mq, which is the largest 
among all synthesized elements. The mass of IMEs (Mg, Si, 
S, and Ca) is ~ O.28M0 in total. The mass of O is esti- 
mated to be O.23M0, and is predominantly located in the 
outer layers. There is no sign of C, and the upper limit is 
M(C) < 0.016M Q at v > 10500 kms" 1 . 



2 These two delayed detonation models have different transition 
densities p = 1 and 2 X 10 7 g cm -3 for WDD1 and WDD2, re- 
spectively. 



From these yields, the validity of the W7 density struc- 
ture can be tested. We use the formula for the kinetic en- 
ergy of SNe la: E K = [1.56M( ° 6 Ni) + 1.74M (stabl e Fe) + 
1.24M(IME) -0.46] x 10 51 erg (|Wooslev et al.ll2007l ). where 
M( 56 Ni), M(stable Fe) and M(IME) are the ejected mass of 
56 Ni, stable Fe-group elements and the IMEs, respectively. 
Substituting the element masses derived for SN 2003du, we 
derive an ejecta kinetic energy of Ek = 1-25 x 10 51 erg, 
which is in a good agreement with the kinetic energy of the 
W7 model (E K ~ 1.3 x 10 51 ). Thus, the element masses de- 
rived from the modelling under the assumption of the W7 
density structure are self-consistent. 

We first compare the integrated yields with those of 
one-dimensional models (Table [5]). The total mass of sta- 
ble Fe-group elements estimated for SN 2003du is reason- 
ably consistent with that of the one-dimensional explosion 
models. The amount of neutron-rich stable elements in the 
central region is sensitive to the deflagration speed, which 
affects the efficiency of the electron captures. Thus, the good 
agreement implies that the deflagration speed in the explo- 
sion models is reasonable. 

The mass of 56 Ni in SN 2003du is also consistent with 
that of the one-dimensional deflagration and delayed deto- 
nation models. In the delayed detonation models, a larger 
amount of 56 Ni suggests an earlier transition to a detona- 
tion. Thus, if the delayed detonation model is the correct 
scenario, a transition density close to that of the models 
WDD1 and WDD2 is reasonable for SN 2003du. 

The total mass of IMEs estimated for SN 2003du is 
smaller than in the delayed detonation models, while the 
mass of O in SN 2003du is larger. This suggests that O 
burning in SN 2003du was not as strong as in the delayed 
detonation models. However, the unburned C mass in the 
outermost layers in SN 2003du is smaller than in the defla- 
gration models. Thus, the outermost part of the C+O WD 
is burned more intensely than in the deflagration model (see 
also Section f6 ■ 1 f) . 

Next, we compare the results with multi-dimensional 
models. The mass of 56 Ni in SN 2003du is clearly larger than 
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Table 2. The total mass of the elements in SN 2003du derived from spectral modelling and nucleosynthesis products of different SN la 
explosion models 



Species SN 2003du W7 1 WDD1 2 WDD2 2 b30.3d.768 3 Y12 4 O-DDT 5 



c 


<1.6E-02* 


4.83E-02 


5.42E-03 


8.99E-03 


2.78E-01 


8.73E-03 


3.49E-03 





2.3E-01 


1.43E-01 


8.82E-02 


6.58E-02 


3.39E-01 


1.07E-01 


1.51E-01 


Na 




6.32E-05 


8.77E-05 


2.61E-05 


8.65E-04 




5.87E-05 


Mg 


2.4E-02 


8.58E-03 


7.69E-03 


4.52E-03 


8.22E-03 


8.70E-2 


2.16E-02 


Si 


2.1E-01 


1.57E-01 


2.74E-01 


2.07E-01 


5.53E-02 


1.27E-1 


2.87E-01 


S 


4.8E-02 


8.70E-02 


1.63E-01 


1.24E-01 


2.74E-02 


7.03E-2 


1.27E-01 


Ca 


1.4E-03 


1.19E-02 


3.10E-02 


2.43E-02 


3.61E-03 


1.82E-2 


1.70E-02 


Ti 




3.43E-04 


1.13E-03 


1.02E-03 


8.98E-05 


1.41E-5 


8.72E-04 


Cr 




8.48E-03 


2.05E-02 


1.70E-02 


3.19E-03 


2.96E-4 


1.04E-02 


Fe a 


1.8E-01 


1.63E-01 


1.08E-01 


1.02E-01 


1.13E-01 


6.50E-3 6 


1.11E-01 


56 Ni 


6.5E-01 


5.86E-01 


5.64E-01 


6.90E-01 


4.18E-01 


9.26E-1 C 


5.40E-01 


Ni d 


2.4E-02 


1.26E-01 


3.82E-02 


5.87E-02 


1.06E-01 




8.05E-02 



*Upper limit at v > 10500 kms . "Not constrained. ***Not given. 
"Stable isotopes except for 56 Fe from 56 Co decay. 
b Only 52 Fe. 

c The sum of Fe-group elements except for 52 Fe, and thus, the upper limit for 56 Ni fsee lPlewalhoOTt) . 

^Stable isotope s ^Ni, 60 Ni, 61 Ni 62 Ni an d 64 Ni 

1 lNomoto etHl dl984r), 2 llwamoto et alj <199Sh. 3 lTravaglio et alj <2004l. 4 |Plewal d2007f). 5 lMaeda et all (|2010TJ| 



in th e three-dimensional deflagration model (|Travaglio et alj 
The three-dimensional deflagration model also pre- 
dicts more unburned C, which is not consistent with the 
small C mass derived for SN 2003du. Thus, also from the 
nucleosynthetic point of view, the three-dimensional defla- 
gration model is unlikely to be a reasonable model for a 
normal SNe la such as SN 2003du. 

It is interesting to compare the yields w ith the deto- 
nating failed deflagration model (|Plewal[2007l) and the off- 
center delayed detonation model ( Maeda et al.ll2010bt ). be- 
cause these models may explain the kinematic offset in the 
central region sugges t ed by the NIR s p ectra at late phases 
(iHoflich et al. | |2004 iMotohara et all l200rj ; iMaeda et al.l 
20100 ). The detonating failed deflagration model (jPlewa 
20071 ) predicts too much 56 Ni compared with SN 2003du. 



Thus, this model does not seem a likely model for nor - 
mal SNe la, as already cautioned by iMeakin et all (|2009f ). 
The yields of the two-dimension al off-center delayed detona- 
tion model (|Maeda et al.|[2010bh are in good agreement with 
those derived for SN 2003du. It is, however, noted that the 
model predicts 56 Ni synthesized near the center, while the 
existence of a region free of 56 Ni is suggested from our one- 
dimensional analysis of the nebular phase spectrum (Section 
0. 



6.3 Comparison with SNe 2002bo and 2004eo 

In this section, we compare the abundance distributions and 
integrated yields derived for SN 2003du with those derived 
for SNe 2002bo (Paper I) and 2004eo (Paper II). SN 2002bo 
is a normally luminous SN (Amis = 1-13 mag), but has high 
line velocities and a high velocity gradient at early phases 
|Benetti et al.ll2004| ). SN 2004eo has a transitional luminos- 
ity between normal and subluminous SNe la, with a rela- 
tively rapid decline o f the light curve (Amis = 1-46 mag, 
iPastorello et alj|2007h . 

These three SNe have in common an innermost layer of 
stable Fe-group elements which are synthesized by electron 



capture. The mass of this region ~ 0. 2Mq is similar in th e 
three SNe studied thus far (see also iMazzali et all 120071 ) . 
although the mass of 56 Ni is different, 0.65, 0.52, and 0.43 
M Q for SNe 2003du, 2002bo, and 2004eo, respectively. 

The most significant difference between SNe 2003du and 
2002bo is in the outermost layers. In SN 2002bo Si reaches 
the outer layers (Paper I). The abundance of Si in the outer 
layers of SN 2002bo is larger than in SN 2003du. This is 
the primary reason why Si line velocities in SN 2002bo are 
hi gher than in SN 200 3du, as demonstrated in Appendix 
[Bl iTanaka et al.l (|200ct ) also suggested that the difference 
between high-velocity and low- velocity SNe la is caused by 
the difference in the Si abundance. 

Alternatively, the difference can also be caused by 
the differenc e in the density structure in the outer layers. 
I Maeda et al.l (|2010al ) pointed that the difference in the ve- 
locity (and/or velocity gradient) is caused by viewing-angle 
effects in an aspherical explosion that has different density 
structures depending on the viewing angle. 

SN 2002bo also has a more extended 56 Ni distribution 
than in SN 2003du. T his explains th e rapid rise of the LC of 
SN 2002bo (Paper IV lPignata et all (|2008l ) suggest that the 
rising time of the LC of SNe with a high velocity gradient 
is generally shorter than that of SNe with a low velocity 
gradient. Thus, the extended 56 Ni distribution may be a 
common property of high-velocity SNe. 

Interestingly, SNe 2003du and 2004eo have a similar 
distribution of Si, O, and C above v = 10000 kms -1 
(M r = 0.9M©). This indicates that the different shapes of 
the LCs of these two SNe result from the difference in the 
inner layers. In fact, a difference between these SNe is seen 
in the mass and mass fraction of 56 Ni. A smaller amount of 
56 Ni is responsible for the fainter luminosity and more rapid 
decline of SN 2004eo (see also Paper II) . 
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7 CONCLUSIONS 

We have studied the element abundance distribution in SN 
2003du by modelling the optical spectra at the photospheric 
and nebular phases. These are used to place constraints on 
the abundances at the outer and inner layers, respectively. 
Since SN 2003du is a normal Type la SN, both photometri- 
cally and spectroscopically, the abundance distribution de- 
rived for SN 2003du can be considered as representative of 
normal Type la SNe. 

The abundance distribution and the integrated yields 
in SN 2003du can be summarized as follows. The innermost 
layers are dominated by stable, neutron-rich elements, with a 
total mass of ~ 0.2Mq. The yield is roughly consistent with 
the results of one-dimensional models. This may imply that 
the deflagration speed assumed in the models is reasonable. 

Above the stable elements, there is a thick 56 Ni-rich 
layer. The total mass of 56 Ni is O.65A/0, which is consis- 
tent with the one-dimensional deflagration model and the 
one-dimensional delayed detonation models with transition 
density ~ 1 — 2 x 10 7 g cm~ 3 . The 56 Ni mass is larger than in 
the three-dimensional deflagration model and smaller than 
in the detonating failed deflagration model. 

The layers dominated by Si and S are located above 
the 56 Ni -rich layers. The dominant element at M r > I.IMq 
(v > f3000 kms -1 ) is O. The progenitor C+O WD is almost 
entirely burned, and the mass of leftover unburned C is quite 
small (< O.O16M at v > f0500 kms" 1 ). The small mass of 
unburned C is consistent with the delayed detonation model. 
However, the mass of O derived for SN 2003du is larger than 
that in the delayed detonation model. These results suggest 
that in the outermost layers C-burning is more intense than 
in the deflagration model, but O burning is not as strong as 
in the delayed detonation models. 

The boundary between stable Fe/ 56 Ni, 56 Ni/IME, 
IME/O-rich layers is not as sharp as in the one-dimensional 
models, suggesting a moderate degree of mixing. However, 
the strong mixing seen in the multi-dimensional deflagration 
model is not preferable. 

This work has made use of the SUSPECT supernova spectral 
archive. V.S. is financially supported by FCT Portugal un- 
der program Ciencia 2008. This research has been supported 
in part by World Premier International Research Center Ini- 
tiative, MEXT, Japan. 
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APPENDIX A: ASSUMPTIONS IN THE CODE 
Al Photospheric Phase 

The spectrum synthesis code for the photospheric phases as- 
sumes a sharply defined spherical photosphere as an inner 
boundary. At the inner boundary, blackbody radiation is as- 
sumed. Since the temperature structure [Tn(r)] is unknown 
at first, the calculation starts with an assumed temperature. 
Using this temperature, the ionization and excitation in the 
optically thin atmosphere above the photosphere are calcu- 
lated considering the dilution of the radiation. The degree 
of dilution W is defined by the equation J = WB(Tr) in 
our code. 

The excitation (the population of an excited level n 7 ; 
j = 1 for th e ground sta te) is computed as in lAbbott fc Lucvl 
(|l985h and lLucvl (|l999h : 



V± = w^e~ tj/kBTR . 
ni gi 



(Al) 



where gj and ej are the statistical weight and the excitation 
energy from the ground level, respectively. 

For the ionization, the nebular approximation, which 
is similar to that made in the studies of stellar winds 
dLucv fc Solomonlll970l; I Abbott fc Lucvlll985h or planetary 
nebulae (jGurz adva nlll997l). are adopted with the modifica- 
tions for SNe (|Mazzali fe Lucvlll993r i: 



N +1 N e 
N 



[S( + (1 - QW]W 



\T R ) V Ni 



(A2) 



where iV e and T e are the electron density and the electron 
temperature, respectively. In the equation, 5 and £ are the 
correction factors for the optically thick continuum shorter 
than the Ca II ionization edge and the fraction of recom- 
binations that go direc tly to the ground state, respectively 
(IMazzali fc Lucvl fl993T) . The first and second terms in the 
first square brackets on the right hand side represents the ex- 
citation from the ground and excited levels, respectively. The 
last term on the right hand side is the ionization computed 
from the Saha equation with Tr. For the electron tempera- 
ture, T e = 0.9Tr is c rudely assumed (|Klein fc Castorlll978l ; 
lAbbott fc Lucvlll985l ). 

With the ionization and excitation states, line opti- 
cal depths are calculated under the Sobolev approximation, 
which is a sound approximation i n rapidly expanding enve - 
lope with large velocity gradients (|Castoilll970l ; lLucvHl97ll ) . 
A number of photon packets are traced above the photo- 
sphere in spherical coordinates, taking into account electron 
scattering and line scattering. Fo r line scattering, the effect 
of the line branching is included (|Lucvll 19991 ; lMazzalill2000l ; 
iPinto fc Eastmanll2000f ). The stream of the photon packets 
gives the flux at each radial mesh and a frequency moment 
is calculated by 



J vj v dv 
J J v dv 



(A3) 



The radiation temperature is estimated from this frequency 
moment via 



hv 
k B T R 



(A4) 



Here x represents for the mean energy of the blackbody ra- 
diation [x = (h/k B T a )( J vB v dv/ J B v dv) = 3.832]. 
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Using this new temperature structure, the ionization 
and excitation are updated. Then, the next Monte Carlo 
ray tracing starts. This procedure is continued until the 
temperature converges. Finally, t he emergen t spectrum is 
calculated using a formal integral |Lucvlll999l ). 

It is known that these assumptions give a good ap- 
proximation to the re sults of detailed NLTE calculations 
l|Pauldrach et al]|l996l ). However, since no energy deposition 
by 56 Ni and 56 Co is assumed above the photosphere, the as- 
sumptions become unreliable as a large part of the 56 Ni-rich 
layer becomes optically thin. In addition, the assumption 
of a sharp photosphere becomes progressively invalid after 
maximum brightness. Thus, the quantities derived by fit- 
ting the spectra at later epochs (especially after maximum 
brightness in SNe la) are more uncertain. 



A2 Nebular Phase 

The spectrum synthesis code for the nebular phases first cal- 
culates the energy deposition by 7-rays and positrons pro- 
duced in the decay of 56 Co. For a given distribution of 56 Ni, 
the energy de position in each shell is evaluated in a Monte 
Carlo scheme (|Cappellaro et alii 19971 ) . The deposited energy 
is spent in the ionization and in the heating of thermal elec- 
trons. The latt er is known to be dominant in Fe-rich ejecta 
(~ 97 per cent. lAxelrodlll980l ). 

For given ionization fractions in the ejecta, the heating 
of the thermal electrons (F) is balanced with the radiative 
cooling (thermal balance): 



r — Tline + Fg 



(A5) 



where Fn nc , Fg and r rcc are the cooling via lines, free-free 
emission, and recombination, respectively. Since the contri- 
bution from the latter two are small er than the line cooling 
(~ 0. 03 and 4 per cent, respectively. iRuiz-Lapuente fc Lucvl 
1 19921 ). we only consider cooling by lines. The line cooling 
per unit volume by lines is given by 



iKotak et all 120051), 2003cg (LVG lElias-Rosa et~ai] 120061 ). 
2005cg (LVG. iQuimbv et all 112000 ). This is in contrast to 
the broad ab s orptio n seen in SNe 20 02bo (Fig. IBTI HVG , 



Benetti et ail 120041). 2002dj (HVG, IPignata et al 



2008 



2006X (HVG. IWang et al.ll2008l ; lYamanaka et al.ll2009l). or a 
possi ble boxy profile i n SNe 1990N (LVG.lLeibundgut et all 
199ll). 2001el /LVG . iMattila et all 12003 ), 2004dt (HVG 



Altavilla et ai]|2007l). 20 05cf (Fig. HU LVG, iGaravini et"ail 
20071 ; I Wang et all |2009| ). Here LVG and HVG mean low 



Time — ^ ^ Tlfcj ^ ^ AjihVij j3ij , 



(A6) 



velocity gradient and high velocity gradient, respectively. 
These are the classifications of Type la SNe based on the 
temp oral evolution of the Si n line velocity jBenetti et alj 
120051 ). Note that high velocity features of Ca 11 lines are 
commonly seen in b oth groups at pre-maximum epochs 
i|Mazzali et al.ll2005al ). 

We use the spectrum at —10.9 days from B maximum 
(fph = 10500 kms" 1 ) to constrain the distribution of Si 
in the outer layers because this spectrum is reproduced 
better than the one at —12.8 days (Fig. [2]|. To reproduce 
the observed profile, a moderately mixed-out Si distribu- 
tion is required. The mass fraction of Si is X(Si)~ 0.3 at 
v — 10500 — 15000 kms" 1 and it is suppressed at the outer 
layers, X(Si)~ 0.02 - 0.05 at v > 15000 kms" 1 (see Figs. [7] 
and|SJ). This model is shown in Figures [3 and [Bl] by the red 
line. 

In the right panel of Figure IBll two other models are 
shown. The orange line ("cutoff") shows a model spectrum 
where the Si abundance has been set to X(Si)=0 at v > 
15000 kms" 1 as in the W7 model (Fig. [8]). In this model, 
the absorption of Si at v > 15000 kms" 1 is clearly too 
weak, suggesting that some degree of burning took place in 
the outer layers (jGerardv et al.|[2004T ). 

The purple line ("homogeneous") shows the model spec- 
trum with a homogeneous Si abundance with A"(Si)=0.3 at 
v > 10500 kms" 1 . In this case, the absorption at a high 
Doppler shift is too strong. This gives a similar profile to 
that in SNe 2002bo and 2005cf. In fact, such an extended Si 
distribution was derived for SN 2002bo in Paper I. 



where njy is the number density of the ion k in the en- 
ergy state j, and Aji and /3ij represents for the Einstein A- 
coefficient and the escape probability of the transition j — > i. 
From this thermal balance, an electron temperature is eval- 
uated in each shell. 

The ionization fractions are determined balancing im- 
pact ionization by non-thermal electrons and recombinations 
(ionization balance). Co llisional ionization by therma l elec- 
trons can be neglected (|Ruiz-Lapuente fc Lucvl ll992l). The 
ionization balance is solved for a given electron temperature. 
Using this new ionization, the thermal balance is reconsid- 
ered, i.e., the temperature and the ionizations are solved 
iteratively. After the convergence, the emergent spectrum is 
computed by integrating the contribution from each shell 
considering the Doppler shift. 



APPENDIX B: SILICON ABUNDANCE AT THE 
OUTER LAYERS 

The Si n line profile in SN 2003du is rather narrow and it is 
similar to that found in, e.g., SNe 2002er (Fig. HU HVG, 



APPENDIX C: UPPER LIMIT FOR THE MASS 
OF UNBURNED CARBON 

Lines of carbon are not clearly seen in the optical spectra of 
SN 2003du. A marginal detection may be se en in correspon- 
dence of the emission pe ak of Si 11 A6355 (|Stanishev et al.l 
120071 : iTanaka et al.ll2008l ). The profile around the emission 
peak is enlarged in Figure IC II and plotted versus Doppler 
velocity measured from the C 11 A6578 line. We use the spec- 
trum at —10.9 days from B maximum (v p ^ — 10500 kms" 1 ) 
to estimate an upper limit for unburned C since the emis- 
sion peak of the Si line is reproduced more nicely than in 
the spectrum at —12.8 days. 

We first divide the ejecta at v > 10500 kms" 1 into 
three parts, with boundaries at v = 15000 and 18000 
kms -1 . At v > 18000 kms" 1 , no constraints can be ob- 
tained on the C abundance because of the low density (weak- 
ness of the line) and the blending with the strong Si 11 line. 
Thus, we could set the upper limit of C as X(C)=0.5 at 
v > 18000 kms" 1 . In practice, since there are other ele- 
ments such as Si and Ca in this velocity space, the mass 
fraction is set to be X(C)=0.4 at v > 18000 kms" 1 in the 
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Figure Bl. The line profile of the Si II A6355 shown in the 
Dopplcr velocity measured from the rest wavelength of the 
Si II line (A = 6355 A). Left: The Si line profile in SN e 
2003du (at —10.9 day s from B maximum , IStanishev et al.ll2007l) . 
2002er (-11.3 day s, iKotak et all |2005|) , 2002b o (-12.0 days 
Benetti et all l2004h . and 2005cf (-11.6 days, iGaravini et al.l 
20071) at similar epochs. These spectra are normalized to SN 
2003du at A = 6355 A. Right: The Si line profile in SN 2003du 
at —10.9 days from B maximum compared with three synthetic 
spectra. Three model spectra are normalized to the best fit model 
at A = 6355 A. The red line shows the best fit model using the 
abundance distribution shown in Figs. [7] and [8] The orange line 
("cutoff") shows the model spectrum where the mass fraction of 
Si, X(Si), is set to be zero above 15000 kms -1 . The purple line 
("homogeneous") shows the model spectrum where the homoge- 
neous Si abundance with Jf(Si)=0.3 is assumed at V > 10500 
kms . 



Figure CI. The line profile around the emission peak of the 
Si II A6355 line shown in the Doppler velocity measured from the 
C II A6578. Left: The profile in SNe 2003du at -10.9 days from 
B maximum compared with the model spectra with X(C)=0.01 
(red, best fit), 0.05 (orange), and 0.4 (purple) at v = 15000-18000 
kms -1 . Right: The profile in SN 2003du at -10.9 days from B 
maximum compared with the model spectra with X(C)=0.002 
(red, best fit) and 0.008 (orange) at v = 10500 - 15000 kms" 1 . 

performed assuming the density structure of W7 model. If a 
delayed detonation mod el was used, X(C) in the outermost 
layers would be smaller (jTanaka et al.ll200ct ). 



best fit model. Since the ejecta mass at v > 18000 kms 1 is 
2.0 x 10~ 2 M Q , the upper limit of the C mass at v > 18000 
kms -1 is 1.0 x 1O~ 2 M . 

The model spectrum with X(C) = 0.05 at v = 15000 - 
18000 kms" 1 shows a noticeable C n line at the emission 
peak of the Si II line (left panel of Fig. IC1[) . Since the ejecta 
mass at v = 15000 - 18000 kms" 1 is 6.0 x 1O" 2 M , the 
upper limit of the C mass in this layer is 3 x 10 _3 Mq, which 
is smaller than the C mass allowed in the outermost layers 
(1.0 x 10" 2 M Q a,t v> 18000 kms" 1 ). 

At v = 10500 - 15000 kms" 1 , the constraint on the 
C abundance is even stronger. The model spectrum with 
X(C) — 0.008 gives a clear absorption trough of the C n 
line (right panel of Fig. [Cl) , Thus, although the ejecta mass 
at v = 10500 - 15000 kms" 1 is large (3.4 x 10 -1 M Q ), the 
C mass at these layers is less than 3 x 10 _3 Mq. 

Thus, most of the C that is allowed must be located 
at v > 18000 kms -1 . Summing up the upper limits to the 
C mass in different layers, we can safely conclude that the 
mass of C ejected in SN 2003du is < 1.6 x 10~ 2 A/ Q at v > 
10500 kms" 1 . This is less than in the deflagration models 
(4.86 x 1O" 2 M ), where X(C)=0.5 at v > 15000 kms" 1 
(see Figs. [7] and [HI and Table [2]). Note that all estimates are 
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